Size-dependent effects on electrical contacts to nanotubes and nanowires 
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Metal-semiconductor contacts play a key role in electronics. Here we show that for quasi-one- 
dimensional (QID) structures such as nanotubes and nanowires, side contact with the metal only 
leads to weak band re-alignement, in contrast to bulk metal-semiconductor contacts. Schottky 
barriers are much reduced compared with the bulk limit, and should facilitate the formation of good 
contacts. However, the conventional strategy of heavily doping the semiconductor to obtain ohmic 
contacts breaks down as the nanowire diameter is reduced. The issue of Fermi level pinning is also 
discussed, and it is demonstrated that the unique density of states of QID structures makes them 
less sensitive to this effect. Our results agree with recent experimental work, and should apply to a 
broad range of QID materials. 
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The early work of Schottky, Mott and Bardeen has laid 
the course for much of the fundamental understanding 
and improvement in the performance of electrical con- 
tacts to bulk semiconductors. However, as new nanoma- 
tcrials are explored for nanoelectronics, the fundamental 
aspects of contacts need to be re-examined due to the 
unique properties of nanostructures. An example is car- 
bon nanotubes (NTs): despite much experimental work, 
it is still unclear whether the contacts are Schottky or 
ohmic, with reports of Schottky contacts for Ti [1] and 
ohmic contacts for Au [2] and Pd [3,4]. However, re- 
cent experimental work [5,6] has suggested that the type 
of contact depends on the NT diameter, with Schottky 
contacts for small diameter NTs and ohmic contacs for 
large diameter NTs. 

From a theoretical perspective, it has been demon- 
strated that Fermi level pinning (crucial in traditional 
semiconductors) is ineffective for quasi-one-dimensional 
nanostructures end-bonded to metals [7]. For NTs side- 
contacted by a metal, modeling has been used to extract 
Schottky barriers from experimental measurements [5], 
but have not addressed the origin of the Schottky bar- 
riers; and atomistic calculations have provided case-by- 
case studies [8-11]. However, a more general theoreti- 
cal understanding for side-contacts to QID structures is 
still missing, especially in light of the recent experimental 
findings. 

In this paper, we present a theoretical and modeling 
analysis of side contacts to nanotubes and nanowires. 
We show that the concepts developed for bulk metal- 
semiconductor contacts do not simply carry over to the 
nanoscale. In particular, band re-alignement due to 
charge transfer is weak due to the limited available de- 
pletion width. In NTs, this leads to relatively small and 
slowly varying Schottly barriers with NT diameter. In 
nanowires (NWs), there is a range of diameters with 
minimized Schottky barriers, providing optimal contact 
properties. We also demonstrate that in general, QID 



structures are much less sensitive to Fermi level pinning 
than their bulk counterparts. Finally, a conventional 
strategy for making ohmic contacts is to heavily dope the 
semiconductor near the contact; we show that at typical 
dopings, the contact resistance increases rapidly as the 
nanowire diameter is decreased. 

We begin by describing the contact geometry consid- 
ered here. Figure la shows a sketch of a cross section of 
the contact consisting of a QID structure embedded in 
a metal. For explicit systems, we consider a single-wall 
NT, as shown in Fig. lb, or a solid nanowire as in Fig. 
lc. For the NT, the metal forms a cylindrical cavity of 
radius R + s where R is the NT radius and s = 0.3 nm is 
the distance between the NT and the metal, while for the 
NW we consider a solid, continuum cylinder embedded 
in a perfect metal, with a sharp interface between the 
nanowire surface and the metal. 

In the simplest picture, the difference between the 
metal Fermi level Ep and the semiconductor valence 
band edge 
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FIG. 1. Panel (a) shows a cross-section of the contact along 
the length of the nanostructure. Panels (b) and (c) show 
radial cross-sections for metal-nanotube and metal-nanowire 
contacts. Panel (d) shows the band alignment before charge 
transfer. In a bulk contact, panel (e), band bending over a 
distance W leads to a Schottky barrier A bu i k . 
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E v (the barrier for holes) is (Fig. Id) 

A = E g + X - $m (1) 

where $ m is the metal workfunction, x is the semicon- 
ductor electron affinity and E g is the semiconductor band 
gap. A positive value for Ao indicates a Schottky bar- 
rier, while a negative value indicates an ohmic contact. 
Because bandgap decreases with increasing diameter for 
Q1D structures, the value of Ao depends on the nanos- 
tructure diameter. The behavior of Eq. (f) for undopcd 
NTs is shown in Fig. 2 for a value of <& m — <&nt — 0.4 eV 
(typical of Pd), and using the relation E g = 2aj/d be- 
tween bandgap and NT diameter d (a — 0.142 nm is the 
C-C bond length, 7 = 2.5 eV is the tight-binding overlap 
integral, and &nt is the NT workfunction assumed to be 
at midgap for an undoped NT). One problem with this 
picture (besides the fact that the physics is incomplete, as 
will be discussed below) is that Eq. (1) predicts large and 
negative values for Ao, signaling strong ohmic contacts. 
However, it is clear that such strong ohmic contacts are 
not observed experimentally. 

In general, charge transfer between the metal and semi- 
conductor leads to band re-alignement. At a bulk semi- 
conductor junction (Fig. le) this charge transfer leads to 
the Schottky barrier 

Abulk =E g + X -$s (2) 

where $ s is the semiconductor workfunction. This rela- 
tionship arises because, in the bulk system, a depletion 
width W perpendicular to the metal-semiconductor in- 
terface is created until the band lineup in Eq. (2) is ob- 
tained. However for Q1D structures, the depletion width 
depends exponentially on the doping [12] and is much 
longer than the device size for non-degenerate doping, 
leading to slowly varying bands outside of the contact; 
and for a three-terminal device the band-bending in the 
channel is governed by the gate voltage. In either case, 
the band alignment 
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FIG. 2. Schottky barrier at nanotube-metal contacts for 
parameters typical of Pd. Dotted line is from Eq. (1), solid 
line is from Eqs (3) and (4) , dashed line is Eq. (6) and circles 
are calculated from an atomistic approach. 



is determined by that in the contact. But for a side- 
contacted Q1D structure, the semiconductor is only a 
few nanometers thick in the direction perpendicular to 
the metal-semiconductor interface; thus only a region of 
the order of the nanostructure cross-section can be de- 
pleted, giving partial band re-alignement. The value of 
A will then be somewhere between Ao and Af, u ik (for an 
undoped NT or NW, At u ik = E g /2, which would always 
give relatively high Schottky barriers). 

Nanotubes are an extreme example of this situation, 
since the possible "depletion width" is the size of the NT 
wall; the charge transfer and image charge in the metal 
create two nested hollow cylinders with opposite charge, 
and an associated electrostatic potential. The charge and 
potential must be self-consistent. We can capture this 
behavior using analytical models for the charge and po- 
tential. The charge per unit area on the NT is expressed 
as 

a = eN J D NT (E + eV NT )f(E - E F )dE (3) 

where Dnt(E) is the NT density of states [13] shifted by 
the electrostatic potential on the NT, f (E — Ep) is the 
Fermi function, and N — 4/(3\/3a 2 ) is the atomic areal 
density. We assume a uniform and sharp distribution of 
the charge on the NT, and all calculations presented in 
this paper are for room temperature. 

For the geometry of Fig. 1, solution of Poisson's equa- 
tion gives the potential on the NT as 

eR R + s e 2 
eV NT = -a — In — — = (4) 

where £0 is the permittivity of free space and C is the ca- 
pacitance per unit area between the metal and the NT. 
Equations (3) and (4) can be solved self-consistently for 
a given NT. In this model the electrostatic potential in- 
duced on the NT modifies the barrier to A = Ao — eV^r- 
Figure 2 shows results of such calculations for parame- 
ters typical of Pd. Clearly, the behavior is different from 
the simple expressions in Eqs (1) and (2). Indeed, the 
bulk limit A}, u ik = E g /2 gives very large barriers, much 
too large to even appear on the scale of Fig. 2. The re- 
sults suggest that there is a transition between Schottky 
and ohmic behavior at a NT diameter around 1.4 nm, 
in agreement with recent experimental data for Pd con- 
tacts [5,6]. We have verified these predictions using an 
atomistic description of the NT based on a self-consistent, 
tight-binding Green's function formalism. As shown in 
Fig. 2, the results of such calculations indicate excel- 
lent agreement with the analytical approach introduced 
above. 

To proceed further we focus on the small and positive 
A regime; approximation of the integral in Eq. 3 leads 
to 
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with (3 = 0.7. Combined with Eq (4) this gives the Schot- 
tky barrier 



(6) 



where a = (e 2 NaV3) / (2y/2j3n 3 / 2 R-(C) . The behav- 
ior of this function is plotted in Fig. 2, showing good 
agreement with the full calculation. The logarithmic de- 
pendence implies relatively slowly varying A, at least 
compared with Eq. (1). The NT diameter delimiting 
Schottky from ohmic behavior is [14] 



d 1 + a 



kT 



NT 



(7) 



where do is the crossover diameter that would be ob- 
tained from Eq. (1). Thus the crossover diameter is 



increased by Sd 



kT 



Making ohmic con- 



tact to a wide range of NT diameters requires a small 
3d; this can be accomplished at low temperature, with 
a large metal workfunction, or with a large capacitance 
(giving a small a). Embedded contacts thus provide an 
advantage over planar contacts because of their larger 
capacitance. 

We now consider side-contacts to nanowires, where the 
possible depletion width increases with diameter, and the 
dependence of the bandgap on diameter is different than 
in NTs. We model a NW with density of states 



D NW {E) = 



V2~r 



(E-EJ2) 



-1/2 



(8) 



where m* is the effective mass. For silicon NWs, it has 
been shown experimentally [15] that the band gap de- 
pends on diameter as E g = Eo + C/d 2 where Eo = 1.12 
eV and C — 4.33 eVnm 2 . We consider the situation 
$m — •ff/if = 0.7 eV typical of contacts to Si. Fig. 3a 
shows the expected Schottky barrier heights from Eq. 
(1) , which predicts ohmic contacts to NWs with diam- 
eters larger than 4 nm. To study the effects of charge 
transfer, we perform a self-consistent calculation of the 
charge and potential, using Eq. (8) to obtain the charge 
and solving Poisson's equation numerically in the NW to 
obtain the potential (we use an atomic volume density 
N v = 5 X 10 28 atoms/m 3 ). 

Fig. 3b, c shows the self-consistent band-bending for 
NWs of 2 and 10 nm radius. Clearly, the nanoscale 
dimension of the NWs prevents the bands from reach- 
ing their asymptotic value; instead, there is only a weak 
band-bending present. To quantify the Schottky barrier 
height, we calculate the spatial average of Ep — E v (r); 
the results plotted in Fig. 3a indicate that the 
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FIG. 3. Panel (a): Schottky barrier at nanowire-metal con- 
tacts for parameters typical of SiNWs. Dotted line is from Eq. 
(1) and solid line is self-consistent calculation. Panels (b) and 
(c): Band-bending across nanowires with diameters of 2 nm 
and 10 nm, respectively; dotted lines are Fermi level. 

contact is always of Schottky character, with the bar- 
rier minimized at a diameter of about 4 nm. Thus, while 
in NTs the barrier height decreases monotonically with 
diameter, the behavior in other Q1D structures may be 
non-monotonic, with a range of diameters providing op- 
timal contact properties. We also note that, just as for 
NTs, the barrier heights are much smaller than the bulk 
limit Ab u ik = Eg/2 (not shown in Fig. 3 for clarity). 

In a bulk metal-semiconductor contact, metal- induced 
gap states (MIGS) lead to Fermi level pinning, and mod- 
ification of the Schottky barrier height to A pin [7]. To 
model this effect in side contacts to Q1D structures, we 
consider a radial pinning charge 



1 pin 



(r) = D N A [Ep - E N (r)} h(r) 



(9) 



where the neutrality level E^ is at midgap [i.e. Ej^(r) = 
-eV(r)], h{r) = e- r / 1 for a NW and h(r) = 5 r>R for 

a NT, and N A = N for a NT and N A = N„ /3 for 
a NW. We choose I — 0.3 nm, a typical value for 
metal-semiconductor interfaces [16]. We add this pin- 
ning charge to Eq. (3) or to the charge calculated from 
Eq. (8) and repeat our self-consistent calculations. 

Figure 4a shows the Schottky barrier calculated for 
several NTs as a function of the density of gap states 
(Api n = Eg/2). Clearly, there is a rapid onset of pin- 
ning at D ~ 0.1 states/ (atom-eV); this value of D is 
rather large considering the van der Waals bonding of 
NTs to metals, and atomistic calculations [8,10] have ob- 
tained seemingly small values. Thus, as in end-bonded 
contacts, we expect that Fermi level pinning will play a 
minor role in side-contacts to NTs. 
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FIG. 4. Schottky barrier as a function of the density of 
gap states for several NTs (a) and NWs (b) . The inset in (b) 
shows the behavior for a planar metal-semiconductor contact. 

Figure 4b shows the effects of Fermi level pinning 
on the barrier height in SiNWs. The results also in- 
dicate a value of D ~ 0.1 states/(atom-eV) required 
to see pinning effects. For comparison, the inset in 
this figure shows the same calculation for a bulk metal- 
semiconductor interface with the same parameters, in- 
dicating that only 0.002 states/ (atom-eV) are needed to 
reach the onset of pinning. Thus, the Q1D system re- 
quires almost two orders of magnitude larger density of 
pinning states compared with the bulk interface. 

The origin of this behavior can be traced to the unique 
density of states of Q1D systems. Indeed, for Si, we 
can repeat the analysis leading to Eq. (5) using the 
density of states for the NW and for the bulk system 
[D bulk {E) = V2(m*f 2 {^n 3 )~ l y/E-E g /2]. This 
leads to the ratio aNwf&buik — (2ttN^ 3 ' f3j / (m*kT). 
The appearance of the kT factor in the denominator is 
entirely due to the Q1D density of states of the NW and 
the presence of a van Hove singularity at the band edge. 
At room temperature, we find that vnw /&bulk > 100; 
thus the MIGS are competing with a much larger charge 
density in the Q1D system. 

Our discussion has so far focused on the situation of 
low doping, where the strategy for making ohmic contacts 
is by selection of a metal with appropriate workfunction. 
In traditional metal-semiconductor contacts, an alterna- 
tive approach is to heavily dope the semiconductor, and 
rely on tunneling through the Schottky barrier to reduce 
the contact resistance and obtain ohmic-like contacts. To 
address the feasibility of this approach for contacts to 
NWs, we repeat our self-consistent calculations for the 
Si NW, focusing on the situation where the metal Fermi 
level is in the middle of the NW bandgap at the inter- 
face, and adding a uniform doping charge of 10 19 e/cm~ 3 . 
Figure 5 shows the band-bending in the presence of this 
doping charge for NWs of 40 and 10 nm diameters. We 
calculate the contact conductance from 

G -Qj^(r%) dE (10) 

where the tunneling probability T(E) is obtained from 
the 
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FIG. 5. Band-bending across Si NWs with doping of 
10 19 e/cm -3 for diameters of 40 nm (a) and 10 nm (b). The 
arrow indicates tunneling of electrons throught the Schottky 
barrier. The normalized resistance is shown in figure (c) as a 
function of NW diameter. 

WKB approximation. The normalized contact resis- 
tance is then Goo/G where Goo is the conductance in the 
limit of large diameters. The behavior of the normal- 
ized resistance as a function of NW diameter is shown in 
Fig. 5c, indicating a rapid increase of the resistance with 
decrease in diameter. The origin of this behavior is the 
increased tunneling distance and reduced range of tun- 
neling energies because of the poor band-bending. One 
implication of this result is that different diameter NWs 
will require different doping levels to achieve the same 
contact quality. 

In summary, we find that the concepts developed to de- 
scribe traditional metal-semiconductor interfaces fail to 
properly account for the properties of contacts to Q1D 
structures. Optimizing device performance will not only 
require selecting Q1D structures for their behavior in the 
channel, but also for their contact properties. We expect 
that our results will be applicable to a broad range of 
Q1D structures. 
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